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INTRODUCTION 


Metal  Matrix  Composites  (MMCs)  are  a  class  of  structural  materials  that  has 
attractive  features  for  many  applications.  Lightweight  MMCs  (based  primarily  on 
the  addition  of  a  reinforcement  to  an  aluminum  or  magnesium  matrix)  are  being 
developed  for  structural  members  such  as  shells,  shelves,  and  stiffeners.  For 
this  generic  type  of  element  the  thickness  is  usually  much  smaller  than  the 
other  dimensions.  Typically,  an  MMC  shell  component  is  no  more  than  2.5mm  thick 
with  the  other  dimensions  being  at  least  one  order  of  magnitude  greater. 

Nondestructive  inspection  of  such  members  frequently  relies  on  ultrasonics 
as  a  tool.  A  geometry  in  which  one  dimension  of  a  part  is  significantly  smaller 
than  the  other  two  is  well  suited  to  the  use  of  plate  modes  of  ultrasonic  wave 
propagation.  With  the  ready  availability  of  ceramic  transducers  that  have 
operating  frequencies  in  the  range  of  1/2  to  1  megahertz  the  inspection  of  thin 
sheet  and  shell  members  by  plate  waves  is  a  practical  approach.  In  addition  to 
flaw  detection,  a  full  set  of  in-plane  material  elastic  constants  may  be 
obtained  from  the  measurements  of  the  proper  plate  wave  velocities.  The  purpose 
of  this  note  is  to  summarize  and  catalogue  the  first-order  equations  which 
express  plate  wave  velocities  as  functions  of  the  elastic  constants  and  density 
of  an  orthotropic  laminate.  In  addition,  the  method  by  which  elastic  constants 
are  determined  from  phase  velocity  measurements  is  reviewed.  In  practice, 
experimental  measurements  often  garner  more  wavespeed  data  than  are  needed  to 
algebraically  solve  the  pertinent  equations.  This  problem  of  over  determination 
of  variables  is  addressed  via  regressive  analysis.  Some  simple  examples  are 
presented  using  unidirectional  boron/aluminum  (B/Al)  and  graphite/aluminum 
(Gr/Al)  as  representative  MMCs. 


REVIEW  OF  THEORY 


It  is  known  that  the  propagation  speed  of  an  ultrasonic  wave  in  anisotropic 
media  is  a  function  of  the  material  elastic  constants  and  density.  For  single 
crystals,  the  elastic  constants  may  be  unambiguously  defined.  For  two-phase 
materials,  such  as  structural  composites,  an  appropriate  model  of  media  response 
to  a  transient  elastic  disturbance  must  be  adopted  in  order  to  quantify  the 
relationship  between  wave  velocities  and  elastic  constants.  When  the  effective 
wavelength  of  the  acoustic  pulse  is  much  larger  than  the  dimensions  of  the 
microstructural  inhomogeneities,  the  effective  modulus  theory  has  been  applied 


for  this  purpose.  Examples  in  the  literature  demonstrate  the  application  of 
effective  modulus  theory  to  relate  ultrasonic  parameters  to  material  properties 
for  a  variety  of  composites.  A  summary  of  the  relevant  equations  for  the  two 
plate  modes  (P  waves  and  SH  waves)  of  present  interest  is  presented  herein. 

The  effective  moaulus  theory  assumes  that  a  composite  material  is  a 
macroscopical ly  homogeneous,  orthotropic  medium.  For  thin  shell  geometries  the 
further  assumption  of  a  plane  stress  state  is  frequently  adopted.  Under  these 
conditions,  the  Hooke's  law  relationship  between  stress  and  strain  is  written 


if  the  stresses  (o],  02.  *12)  and  the  strains  (ci,  e2»  y12)  are  indexed  to  a 
cartesian  coordinate  system  aligned  with  the  laminate  principal  axes. 

The  problem  of  determining  the  velocity  of  waves  propagating  in  an 
isotropic  plate  was  originally  solved  by  Lamb.  Among  other  things,  the  results 
showed  that  two  basic  families  of  waves  could  propagate,  one  having  symmetrical 
and  one  unsymmetrical  characteristics.  The  long-wavelength  (where  the  wave¬ 
length  is  measured  in  units  of  plate  thickness)  symmetrical  mode  is  mildly 
dispersive,  with  the  approximate  wavespeed  being  a  simple  function  of  the  plate 
moduli  and  density.  The  antisymmetrical  mode  is  highly  dispersive.  For  the 
symmetrical  Lamb  wave,  as  the  wavelength  becomes  infinitely  long,  the  limiting 
velocity  is  referred  to  as  the  extensional  plate  velocity.  Subsequent  work  by 
other  investigators  developed  similar  plate  velocity  expressions  for  orthotropic 
elastic  plates.  The  extensional  plate  wave  velocity  and  its  shear  mode 
analogue  may  be  expressed  as  functions  of  the  elastic  constants  of  a  composite 
laminate  by  using  Equation  (1)  in  conjunction  with  the  equations  of  motion  of 
the  laminate  material  as  it  is  disturbed  by  the  ultrasonic  wavefront.  A  simple 
derivation,  based  on  a  time-harmonic  wave  train,  is  summarized  here. 

A  field  point  in  the  laminate  located  at  coordinates  (x,y)  oscillates  with 
respective  displacements  (u,v).  The  strains  at  (x,y)  are  then 
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When  the  wavelength  of  the  disturbance  is  much  greater  than  the  laminate 
thickness  plate  wave  modes  are  activated  and,  for  the  lowest  order  of  P  wave  and 
SH  wave,  the  displacement  at  (x,y)  may  be  approximated  as  uniform  across  the 
plate  thickness.  The  displacements  due  to  a  harmonic  wave  will  vary  with  time, 
t,  as 


u  *  u0  exp[i(k  x  cose  +  k  y  sine  -  wtj 


(3) 


A  similar  expression  exists  for  v.  In  Equation  (3),  u  is  the  angular  frequency 
k  *  u/V,  where  V  is  the  phase  velocity  of  the  wave,  and  a  is  the  angle  between 
the  x  axis  and  the  normal  to  the  wavefront. 

From  Newton's  laws,  the  equation  of  motion  of  the  field  point  is 


In  Equation  (4),  p  is  the  composite  density. 

Combining  Equations  (1)  through  (4)  yields  the  Christoffel  equation  for 
(u0,  vo)»  which  is  solved  by  setting  the  characteristic  determinant  equal  to 
zero.  This  relationship  is  known  as  the  dispersion  equation.  Equation  (5). 


A11  “  X 


A22  -  X 


In  Equation  (5),  X  *  pV2  and  the  A-jj  parameters  are  functions  of  the  laminate 
stiffness  coefficients,  Q<,-,  and  the  direction  of  wave  propagation,  e,  as 
follows: 

An  ~  Qll  cos2e  +  Q66  sin2e 

a22  =  Q66  cos2e  +  Q22  sin2e  ( 

A12  s  (Q12  +  Q66)  sine  cose 

Mathematically,  Equation  (5)  is  a  quadratic  with  two  roots  for  each  angle 
e.  The  larger  root,  representing  the  faster  wave,  Xp,  is 


2  +  4  A2 


The  smaller  root,  corresponding  to  the  slower  wave,  Xs,  is  obtained  by  changing 
the  plus  sign  between  the  two  terms  of  Equation  (7)  to  a  minus  sign. 

The  physical  interpretation  of  the  two  roots  is  that  in  any  direction  in 
the  plane  of  the  lamina  two  types  of  waves  may  propagate  which  involve 
principally  in-plane  stresses.  Neither  type  will  generate  a  net  moment  as,  for 
example,  does  a  wave  of  flexure.  These  waves  become  pure  mode  longitudinal  and 
pure  mode  shear  only  when  the  propagation  direction  is  along  a  laminate 
principal  axis.  For  these  directions  e  *  0°  or  e  *  90°  and  the  phase 
velocities  are 
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VF(0°)  *  /tln/p 

Vp(90°)  *  /q22/p  (8) 

Vs(0°)  •  Vs(90°)  =  ^66/p 

For  off-axis  waves  the  mode  is  not  pure  and  the  terms  pseudo-longitudinal  waves 
and  pseudo-shear  waves  are  sometimes  used  as  descriptors.  The  relationships 
between  the  Qij  coefficients  and  the  conventional  engineering  elastic  constants 
(Young's  moduli,  shear  modulus,  and  Poisson's  ratios)  are: 

Qll  =  Ell/(1  -  vi2  v2i) 

Q22  =  E22/(l  -  vi2  v2l) 

Ql2  =  v]2  Q22  (9) 

Q66  ■  g12 

v21  En  =  vi2  E22 

In  Equation  (9),  En  and  E22  are  the  extensional  moduli  in  directions  x 
and  y,  respectively;  vi2  and  v?i  are  the  major  and  minor  Poisson's  ratios;  and 
G12  is  the  in-plane  shear  modulus. 

In  order  to  check  the  predictions  of  these  equations  a  panel  of  orthotropic 
material  with  known  elastic  constants  is  needed.  Given  the  set  of  four  elastic 
parameters  plus  the  panel  density,  the  phase  velocities  of  the  two  types  of 
waves  may  be  calculated  from  Equation  (8)  and  compared  to  the  experimental 
measurements  in  the  principal  directions.  Off-axis  waves  may  also  be  used  for 
this  purpose  but  greater  care  must  be  taken  to  ensure  that  the  correct  mode  is 
observed  due  to  certain  directional  propagation  effects  present  in  orthotropic 
media,  as  discussed  below. 

Acoustical  birefringence  occurs  in  anisotropic  media.  That  is,  the  phase 
velocity  varies  with  direction.  Even  when  using  a  model  of  material  behavior 
that  neglects  dependence  of  phase  velocity  on  frequency,  the  directional 
variation  of  V  leads  to  the  existence  of  a  wave  modulation  envelope  which 
travels  with  group  velocity,  Vq,  that  is  not  always  the  same  as  the  phase 
velocity.  The  group  velocity  is  a  vector  and  its  components  in  the  x  and  y 
directions  are 

(10) 
"  S7 

In  Equation  (10),  kx  =  k  cose  and  k y  =  k  sine.  The  partial  derivatives 
are  determined  by  implicit  differentiation  of  the  dispersion  relation.  Equation 
(5).  This  relation  written  in  expanded  form  is 
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o(w,kx,ky)  =  (kx  2Qi i  +  ky  2Q66  -  pw2)  (kx  2Q66  +  ky  2Q22  -  pw2) 

-  (kx  ky(Qi2  +  Q66))2  =  0 
The  components  of  Vg  may  be  expressed 

Vgx=-^ 

bX  90/aw 


Vgy  =  - 


ao/ak^ 
30/  3a» 


The  corresponding  partial  derivatives  are 

s  (kx  2Qll  +  ky  2Q66  -  P«2) (2  kx  Q66) 

+  (kx  2Q66  +  ky  2Q22  -  pw2)(2kx  Qn)  -  2kx  ky  2(Qi2  +  Q66)2 


=  (k,  2Qn  +  ky  2q66  -  p(i»2)(2  ky  Q22) 


+  (kx  2Q66  +  ky  2Q22  -  pw2)(2ky  Q66>  -  2kx  2  ky  (Q12  +  Q66)2 


—  *  (kx  2(Qll  +  Q66>  +  ky  2 (Q22  +  Q66)  -  2pw2)(-2p»)  (14) 

3w 

The  direction  of  group  velocity  advancement  with  respect  to  the  x  axis,  *g, 
is  found  from 


The  magnitude  of  Vg  is 


(16) 


In  Equation  (15),  eg  is  the  angle  of  inclination  of  the  group  velocity 
propagation  direction  for  an  ultrasonic  wave  packet  that  is  launched  with  a 
phase  velocity  propagation  direction  of  angle  e.  Solution  of  Equation  (15) 
shows  that  eg  *  e  only  for  the  two  angles  0°  and  90°. 

The  fact  that  for  off-axis  waves  the  phase  and  group  velocities  do  not 
advance  In  the  same  direction  is  important  to  properly  sizing  a  specimen  for 
measurement  of  off-axis  phase  velocities.  When  a  pulse  of  acoustical  energy  is 
radiated  by  a  plane  wave  transducer,  the  phase  front  is  aligned  parallel  to  the 
transducer  face.  However,  the  energy  propagates  in  the  direction  of  the  group 


Vg  *  [Equation  (12)]2  +  [Equation  (13] 
[Equation  (14)]^ 
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velocity,  which  is  not  normal  to  the  transducer  face.  The  angle  of  energy  cant 
from  the  transducer  normal  may  be  significant  and,  for  a  through-transmission 
type  of  wave  speed  measurement,  can  lead  to  improper  placement  of  the  receiving 
transaucer.  The  specimen  dimensions  must  be  large  enough  to  account  for  this 
canting  effect  so  that  the  receiver  may  be  offset  by  the  proper  amount  to 
intercept  the  pulse.  Figure  1  depicts  schematically  the  beam  deflection 
phenomenon  for  an  off-axis  wave  propagating  in  a  directional  media.  From  the 
figure,  the  phase  velocity  is  numerically  equal  to  the  perpendicular  distance 
separating  the  transducers,  L,  divided  by  the  time  of  flight  of  the  group  wave 
packet.  Also,  note  that 

V  =  VG  cos(eG  -  e)  (17) 

A  number  of  investigators  have  used  measured  phase  velocities  to  calculate 
the  elastic  constants  of  a  composite.  In  the  present  formulation,  this  problem 
reduces  to  determining  the  four  stiffness  parameters--Qn ,  Q22.  Q66»  Ql2— from 
measured  plate  velocities.  Therefore,  a  minimum  of  four  velocities  must  be 

recorded.  At  least  one  value  is  needed  in  an  off-axis  direction  or  vl2 
cannot  be  determined.  Algebraic  solution  of  Equation  (7)  from  the  four 
experimental  quantities  plus  the  density  is  then  possible.  A  convenient 
approach  is  to  interrogate  a  specimen  at  e=0°  and  e=90°,  since  for  these 
directions  the  phase  velocity  propagates  normal  to  the  transducer  face.  The 
influence  of  Poisson's  ratio  on  wave  speed  derives  largely  from  the  effect  of 
the  magnitude  of  A12  on  the  calculated  results  of  Eq^tion  (7).  The  value  of 
A]2  is  proportional  to  the  product  sine  cose,  which  '  a  maximum  at  e=45°. 
Therefore,  it  is  advisable  to  determine  at  least  one  wave  speed  at  e=  */4, 
if  possible. 

At  times,  more  than  four  experimental  wave  speed  values  are  available. 

When  this  situation  occurs,  the  accuracy  of  the  Q -j j  determinations  may  be 
increased  by  using  the  method  of  non-linear  least  square  error  analysis.  The 
approach  is  to  fit  the  two  curves  defined  by  Equation  (7)  to  the  total  set  of 
data  in  such  a  way  that  the  sum  of  the  absolute  deviations  of  the  data  points 
from  the  smooth  curves  is  a  minimal.  The  following  discussion  assumes  that  the 
experimental  error  in  wave  speed  measurement  is  the  same  for  all  data  points. 

From  Equation  (7),  two  curves  represent  the  behavior  of  both  in-plane  plate 
wave  modes,  where  the  curve  parameters  are  Xf  and  Xs  expressed  as  functions 
of  the  four  Q-jj's  and  e.  The  goodness  of  fit  of  a  set  of  data  to  these  two 
curves,  «,  is  defined  as 

6  =  ^  {(XF  -  X1>*  +  <xs  -  x2)2}  08) 

Equation  (18)  is  a  summation  over  all  specimen  angles  interrogated,  n.  In 
Equation  (18),  Xp  =  Xp(Qij,en)  and  X]  =  pV^,  where  Vp  is  the  measured  speed  of 
the  faster  wave  for  angle  en.  Similarly,  Xs  =  Xs(Q-jj,en)  and  X2  =  pV §  at  angle 
®n* 
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The  function  a  is  here  considered  a  continuous  function  of  the  four 
unknown  Qij's;  each  Qij  is  an  independent  variable.  Thus  a  describes  a  hyper¬ 
surface  in  four-dimensional  Q^j  space  and  the  space  must  be  searched  for  the 
minimum  value  of  a.  The  searcn  procedure  consists  of  establishing  a  four¬ 
dimensional  grid  mapping  of  the  area  of  Qij  space  in  which  a  is  suspected  to 
be  minimal,  and  then  progressively  tightening  the  grid  by  reducing  the  range  of 
values  of  each  of  the  Qij's. 

In  the  simplest  mapping  procedure,  the  permissable  range  for  each  Qij  is 
divided  into  two  equal  increments  about  its  median  value.  The  mid-range  point 
plus  the  two  extremities  constitute  three  trial  values  for  each  Qij  so  that  the 
4-parameter  space  is  divided  into  a  grid  consisting  of  81  vertices.  The  value 
of  a  is  then  calculated  at  each  of  these  vertices.  This  procedure  yields  a 
map  of  the  behavior  of  a  as  a  function  of  all  of  the  Qij  parameters.  The  grid 
may  be  made  initially  very  coarse  on  a  "best  guess"  basis,  and  then,  when  a 
local  minimum  is  identified,  made  finer  by  adjusting  the  median  value  of  each 
Qij  and  reducing  the  incremental  offset  about  the  mid-point. 

Computer  programs  in  BASIC  for  calculating  the  phase  and  group  velocities 
from  laminate  elastic  constants  and  density  and  for  determining  the  best-fit  set 
of  Qij's  from  a  collection  of  experimental  phase  velocities  are  included  in 
this  report  as  Appendices  A  and  B,  respectively. 


EXPERIMENTAL 


Two  types  of  MMC  were  available  for  preliminary  evaluation  of  the 
preceeding  theory.  Both  consisted  of  undirectionally  reinforced  aluminum  in  a 
thin  lamina  form  and  are  therefore  characterized  with  respect  to  in-plane 
stiffness  behavior  by  four  independent  elastic  constants— E] ],  the  axial  elastic 
modulus;  E?2»  the  transverse  elastic  modulus;  vj2»  the  major  Poisson's  ratio; 
and  G]2»  the  in-plane  shear  modulus.  Some  physical  properties  of  the  as-received 
plates  are  given  in  Table  1.  The  boron  filament  used  in  B/Al  is  a  product  of 
AVCO;  the  precursor  wire  of  which  the  plate  of  Gr/Al  is  made  was  manufactured  by 
Material  Concepts,  Inc.  Plate  consolidation  was  performed  by  DWA  Composite 
Specialties  in  both  cases.  Although  similar  in  regard  to  matrix  type  and 
reinforcement  arrangement,  the  boron  filament  and  graphite  fiber  differ  in 
average  diameter  by  approximately  one  order  of  magnitude.  Thin  panels  of  these 
kinds  can  serve  as  structural  materials  in  their  undirectional  form  or  as  the 
unit  of  construction  for  complex  laminate  layups. 

The  experimental  approach  was  to  cut  the  panels  into  various  sample  sizes 
for  measurement  of  ultrasonic  phase  velocities  by  through-transmission  methods. 
Figure  2  shows  three  typical  samples  of  B/Al.  Both  longitudinal-type  waves  and 
in-plane  shear  waves  were  studied  using  broadband  transducers  which  had  a 
nominal  0.5  MHz  center  frequency.  The  materials  were  interrogated  to  determine 
the  velocities  of  both  types  of  waves  at  five  angles  with  respect  to  the 
reinforcement  direction.  Densities  were  obtained  by  immersion.  The  resulting 
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TABLE  1.  DESCRIPTION  OF  PLATES  USED  IN  ULTRASONIC  STUDIES 


Boron/Ala 

Graphite/Alb 


Boron/Ala 

Graphite/Alb 


Length 

Width 

Thickness 

(cm) 

(cm) 

(cm) 

23.8 

7.6 

0.138 

30.0 

30.0 

0.068 

Filament 

Alloy 

Reinforcement 

Density 

Type 

Fraction 

(gm/cmr) 

5.6-mil  boron 

6061 

0.51 

2.59 

GY-70  graphite 

A201 

0.37 

2.52 

a.  DWA  Plate  No.  B-1379-1 

b.  DWA  Plate  No.  G-5100 
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ARROW  DENOTES  NORMAL  DIRECTION  TO  TRANSDUCER  FACE 


FIGURE  2.  BORON/ALUMINUM  SPECIMENS 
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velocities  and  calculated  values  of  Xp  and  Xs  are  presented  in  Tables  2  and 
3.  The  composite  elastic  constants  were  obtained  for  each  material  by  applying 
the  least-square-error  program  of  Appendix  B.  These  values  are  listed  in  Tables 
4  and  5. 

The  B/Al  plate  was  not  large  enough  to  permit  additional  mechanical  tests 
to  be  conducted  which  would  provide  comparative  data  for  the  elastic  constants. 
Therefore,  the  results  of  an  investigation  in  the  literature  is  utilized  for 
this  purpose,  Reference  1.  In  this  study  the  complete  orthotropic  elastic¬ 
stiffness  matrix  of  undirectional ,  Borsic-fi lament-reinforced  A1  composites  was 
experimentally  evaluated  for  three  different  volume  fractions  by  ultrasonic  bulk 
wave  velocity  measurements.  The  fiber-volume  fractions  reported  are  0.13,  0.34 
and  0.54.  Interpolation  of  the  elastic  constants  from  this  paper  to  the  fiber- 
fraction  of  the  B/Al  plate  used  in  the  present  study  (0.51)  gives  the  set  of 
elastic  constants  listed  in  Table  4. 

The  Gr/Al  plate  was  large  enough  to  allow  coupons  to  be  cut  for  tensile 
evaluation  of  Young's  modulus  in  the  fiber  direction,  E]].  Three  tests  were 
performed;  the  results  are  given  in  Table  5.  Insufficient  material  precluded 
running  mechanical  tests  to  determine  the  remaining  three  properties.  Since 
£22  and  Gl2  are  both  matrix-dominated  properties,  a  rough  estimate  of  their 
value  for  an  aluminum  composite  incorporating  37 %  graphite  may  be  obtained  from 
studies  on  other  types  of  unidirectional  Gr/Al  of  approximately  equal  fiber 
fraction  (References  2  and  3). 

The  composite  cylinder  assemblage  model  has  been  used  to  predict  the 
properties  of  a  Gr/Al  material  from  calculations  based  on  fiber  and  matrix 
constituent  properties.  Reference  2.  Fiber  fractions  ranging  from  0  to  50 
percent  were  addressed  in  10  percent  increments,  allowing  interpolation  to  the 
value  of  0.37.  Ultrasonic  evaluations  based  on  bulk  waves  of  a  30  volume 
percent  Gr/Al  composite  were  also  reported.  Reference  3.  Both  studies  involved 
graphite  filaments  of  another  type  than  GY-70.  The  elastic  constants  from 
these  papers  are  reproduced  in  Table  5. 

Predictions  of  ultrasonic  group  velocity  by  the  effective  modulus  theory 
were  in  part  evaluated  by  an  experiment  which  delineated  the  group  propagation 
envelope  directly.  A  small  source  of  extensional  waves  was  applied  to  one  edge 
of  the  GY70/A201  plate,  as  shown  in  Figure  3.  A  second  small  transducer  was 
employed  to  track  a  constant  phase  point  in  polar  fashion  about  the  disturbance 
source  as  origin.  This  procedure  yields  a  plot  of  V§  versus  og  for  angles 
ranging  from  0°  to  85°.  The  results,  presented  in  Figure  4,  correspond  to  the 
locus  of  the  wavefront  for  the  faster  plate  wave  at  a  time  of  t  =  12.8  nsec 
after  the  pulse  was  initiated.  The  theoretical  curve,  calculated  from  Equation 
(15)  and  Equation  (16)  using  the  GY70/A201  moduli  of  Table  5,  is  presented  on 
the  same  figure  as  the  experimental  data  for  comparison. 


*The  filament  axial  elastic  modulus  was  approximately  50  x  10^  psi  in  both 
studies  cited.  The  axial  modulus  of  GY70  is  approximately  70  x  10°  psi. 
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TABLE  2.  BORON/ ALUMINUM  UNIDIRECTIONAL  LAMINA  PLATE  WAVE 
VELOCITIES  AND  ACOUSTICAL  PRESSURES 


e 

Vp(mm/psec) 

Vs(mm/ysec) 

XF(GPa) 

Xs(GPa) 

0° 

9.624 

4.741 

239.52 

58.126 

20° 

9.207 

4.859 

219.21 

61.055 

45° 

8.142 

5.300 

171.43 

72.641 

60° 

7.946 

5.274 

163.28 

71.930 

90° 

7.826 

4.820 

158.38 

60.079 

TABLE  3.  GRAPHITE/ALUMINUM  UNIDIRECTIONAL  LAMINA  PLATE  WAVE 
VELOCITIES  AND  ACOUSTICAL  PRESSURES 


e 

VjJmm/iisec) 

Vs(mro/usec) 

XF(GPa) 

Xs(GPa) 

0° 

8.84 

2.74 

196.9 

18.92 

30® 

7.49 

2.90 

141.4 

21.19 

45® 

6.41 

2.96 

103.5 

22.08 

60® 

5.06 

3.11 

64.52 

24.37 

90® 

3.88 

2.75 

37.93 

19.06 
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TABLE  4.  BORON/ ALUMINUM  UNIDIRECTIONAL  LAMINA  ELASTIC  CONSTANTS 


Plate  Waves4  Ref.  lb 


Ell(msi) 

32.8 

32.6 

E22(msi) 

21.9 

19.3 

6i2(msiJ 

8.2 

8.4 

v12 

0.24 

0.25 

(a) 

51  v/o  Boron/ 6061 

(b) 

52  v/o  Bor si c/6061 

TABLE  5. 

GRAPHITE/ALUMINUM 

UNIDIRECTIONAL 

LAMINA  ELASTIC  CONSTANTS 

Plate  Waves4 

Tensile  Test4 

Ref.  2b 

Ref.  6C 

L] i(msi ) 

27.4 

27.0,  29.3,  27.6 

26.6 

23.3 

E22 (rosi ) 

5.0 

- — 

5.1 

4.3 

Gl2(msiJ 

2.5 

— •» 

3.1 

2.8 

U.26 

— 

0.34 

0.28 

(aj  37  v/o  GY70/A201 

(b)  mathematial  model  predictions  for  37  v/o  P55/6061 

(c)  30  v/o  T50/A201 
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LENGTH  (cm) 

FIGURE  4  GROUP  VELOCITY  ENVELOPE  FOR  WAVES  IN  GY70/A201 


DISCUSSION 


A  comparison  of  the  elastic  properties  of  B/Al  and  Gr/Al  from  ultrasonic 
plate  waves  with  values  found  in  the  literature  shows  reasonable  agreement  for 
the  majority  of  the  constants.  The  discrepancies  can  be  attributed  to  sample 
differences  between  this  study  and  the  references  cited  and  the  effect  of 
dispersion  on  plate  wave  velocity  determinations.  The  good  agreement  of  plate 
wave  values  of  Eh  with  the  results  of  other  types  of  measurement  is 
particularly  encouraging,  in  that  axial  performance  of  a  unidirectional 
composite  is  usually  of  first  order  importance  for  design  considerations. 

Reasonable  agreement  is  also  observed  between  the  predicted  and  observed 
group  velocity  envelope  of  an  extensional-type  wave  generated  by  a  point  source 
of  disturbance  in  Gr/Al.  This  fact  lends  confidence  to  the  hypothesis  that 
effective  modulus  theory  can  be  applied  to  relate  the  material  stiffness 
parameters  of  a  unidirectional  MMC  to  elementary  elastic  wave  analysis. 


CONCLUSIONS 


Effective  modulus  theory  appears  adequate  to  provide  a  first-order 
approximation  of  in-plane  plate  wave  propagation  characteristics  in  thin  panels 
and  shells  of  unidirectional  B/Al  and  Gr/Al.  The  plate  wave  approach  appears 
especially  relevant  to  determining  the  value  of  the  major  elastic  modulus. 

Eh,  which  is  of  much  interest  to  both  users  and  manufacturers  of  MMCs. 


RECOMMENDATIONS 


Additional  plate  wave  studies  of  the  present  type  should  be  performed  on 
other  types  of  MMCs  and  on  layup  arrangements  other  than  unidirectional  plies. 
The  effect  of  residual  stresses  resulting  from  the  manufacturing  process  on 
phase  velocities  needs  to  be  investigated.  The  development  of  a  model  of  MMC 
material  behavior  with  regard  to  transient  elastic  disturbances  that 
incorporates  dispersion  should  be  addressed. 


i 
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APPENDIX  A 

COMPUTER  PROGRAM  FOR  CALCULATING  PLATE  WAVE  SPEEDS 
FROM  ELASTIC  CONSTANTS 


The  following  computer  program,  written  in  BASIC,  calculates  the  fundamental 
symmetric  Lamb  and  SH  shear  plate  wave  velocities  of  an  orthotropic  lamina  from 
its  four  elastic  constants  and  density.  Enter  the  input  data  in  Step  2000  in 
the  order: 


ill  (Msi);  E22  (Msi);  G]?  (Msi);  v]2J  and  p  (gm/cm3).  A  sample 
nation  is  presented  in  Table  A-l  using  En  s  27.4  Msi;  Epp  *  5. 


E 

calculation  is  presented 
G ip  *  2.5  Msi;  vi2  *  0.26;  and  p  =  2.52  gm/cm' 
set  of  values  for  GY  70/A201  Gr/Al. 


i1 


E22  *  5.0  Msi; 
which  is  a  representative 


A-l 
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10  REM  PROGRAM  TO  CALCULATE  UAVE  SPEEDS  OF  FUNDAMENTAL  LAMB  MODE 
12  REM  AND  SH  SHEAR  NODE  AS  A  FUNCTION  OF  AN6LE  FOR  THIN  PLATE  MATERIAL 

14  REM  USING  THE  MODULI  AS  INPUTS.  INPUT  Ell ( MSI ) v E22< MSI > v 61 2< MSI > 

16  REN  POISSONS  RATIO,  AND  DENSITY  (6/CC) 

17  PRINT 

15  PRINT 

20  PRINT  "PR06RAM  CALCULATIONS  MADE  AT  S  DE6REE  INTERVALS" 

AO  DIM  T3(95),X1 (95),X2(95),Y1 (95),Y2(95) 

AS  DIN  T9(95),N9(95) ,22(95) ,Z3( 95) 

90  PRINT 

102  READ  E1,E2,6,V1,R 
105  PRINT 

10A  PRINT  ”  Ell  *  ",E1,"MSI" 

108  PRINT  "  E22  «  ",E2,*MSI" 

110  PRINT  "  612  «  •,S,"MSI" 

112  PRINT  "  POISSONS  RATIO  *",V1 

114  PRINT  "  DENSITY  *  ",R,"6HS/CCM 

115  PRINT 
1 1 A  PRINT 

117  PRINT  “FAST  UAVE  SPEEDS  VERSUS  ANGLE" 

118  PRINT 

119  PRINT 

120  PRINT  "THETA  PHASE", "F  PHASE  VEL", "THETA  GROUP", "F  GROUP  VEL" 

122  PRINT  "(DEGREES)", "<NH/USEC)", "(DEGREES)", "(MM/USEC)" 

124  PRINT 

128  V2*V1*(E2/E1 ) 

130  A1*E1/( ( 1 -V1*V2)*0. 1 45) 

140  A2*E2/((1-V1*V2)*0.145) 

150  A3=(V2*E1 )/((1-V1*V2)»0.145) 

UO  A6*G/.145 

200  FOR  1*0  TO  90  STEP  5 

210  T2*I*( 6.2831 85/360) 

215  T3(I)*T2 
220  L*C0S(T2) 

230  L2*L*L 

240  N»(1-L2)**0.5 

250  M2«M*M 

260  H1«L2*A14H2*A6 

270  H2«L2*A64H2*A2 

280  H3*L*H*(A3+A6) 

290  H4*H3*H3 
300  H7«H1*H2 
310  H8*H7*H7 
320  H9*H1*H2-M4 

400  X1-0.5*(H7*(H8-4*H9)**0.5) 

410  X2*0.5*(N7-(H8-4*H9)*40.5) 

420  U1«(X1/R)M0.5 

AOO  F 1 ■ ( L2* A 1 «M2*A6-R*U 1 *U 1 ) • ( 2*M* A2 ) 

A10  G1-(L2*A14M2*A6-R*U1*U1 )*(2*L*A6) 

A20  F2* (L2M6*M2*A2-R*U1 *U1 ) * ( 2*M*A6 ) 

A30  62* ( L2*A6+H2*A2'R*U 1 *U 1 ) • ( 2*L*A1 ) 

AOO  F3«-2*H*L2*(A3*A6)**2 
A50  83«-2*L«M2*(A34A6)**2 
AAO  F4-FHF24F3 

A-2 


670  64*61+62+63 

610  U1*L2*(A1 +A6)+H2*(A6+A2)-2*R*U1*U1 
690  U2*-2*R*U1 
700  II3*U1*U2 

710  M4* ( F4+F4+G4+64 ) / ( U3*U3 ) 

720  U5*U4**.S 
730  T7*F4/64 

740  T8* ( 180/3. 141 59)*ATN(T7) 

760  PRINT  I,U1,T8,U5 
1800  NEXT  I 
1010  PRINT 
1020  PRINT 

1022  PRINT  "SLOW  HAVE  SPEEDS  VERSUS  ANGLE" 

1024  PRINT 
1026  PRINT 

1030  PRINT  "THETA  PHASE ","8  PHASE  VEL", "THETA  GROUP" f"S  6R0UP  VEL11 
1032  PRINT  " ( DEGREES ) " f " ( NH/USEC ) " , " ( DE6REES )  * , "  ( NH/USEC ) " 

1034  PRINT 

1200  FOR  1*0  TO  90  STEP  S 
1210  T2«I*<6. 283185/360) 

1215  T3(I)*T2 
1220  L  »C0S(T2) 

1230  L2«L*L 

1240  H*(1-L2)*+.5 

1250  H2*N*H 

1260  H1-L2+A1+H2+AA 

1270  H2*L2«A6+N2*A2 

1280  H3*L*H*(A3+A6> 

1290  H4=H3*H3 

1300  H7*H1+H2 

1310  H8=H7*H7 

1320  H9*H1*H2-H4 

1400  X1«0.5*(H7+<H8-4*H9)«*.5> 

1410  X2*0.5+(H7-(H8-4*H9)**.5) 

1430  U1*(X2/R)**.5 

1600  F1*(L2*A1+H2*A4-R*U1*U1><<2*N*A2> 

161 0  61  * (L2+A1 +H2+A6-R+U1 *U1 >•( 2+L+A6 ) 

1620  F2*(L2*A6+H2*A2-R«U1«U1 )*<2*H*A6) 

1630  62«<L2*A4+H2*A2-R+U1*U1)*<2+L*A1> 

1640  F3*-2+H«L2*(A3+A6)**2 
1650  63*’2*L*H2*(A3+A6)**2 
1660  F4*F1+F2+F3 
1670  64*61+62+63 

1680  HI *L2* (A1+A6) +H2* ( A6+A2 ) -2*R*U1 *U1 
1690  U2«-2*R*U1 
1700  U3*U1*U2 

1710  U4*(F4*F4+64*64)/(U3*U3) 

1720  U5*U4**.S 
1730  T7-F4/64 

1740  T8*(180/3.14159)MTN(T7> 

1770  PRINT  I,U1,T8,U5 
1780  NEXT  I 

2000  BATA  27.4,5,2.5,. 26,2.52 
9999  ENB 


,  »  •  * 
VWV 
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TABLE  A-l.  CALCULATED  PLATE 

WAVE  SPEEDS  OF  GY70/A201 

Ell  «  27.4 

HS! 

E22  *  5 

RSI 

S12  *  2.S 

HSI 

POISSONS  RATIO’ * 

.24 

OENSITT  ■  2.52 

SHS/CC 

FAST  HAVE  SPEEDS  VERSUS  ANGLE 

THETA  PHASE  F  PHASE  VEL 

THETA  GROUP 

F  GROUP  VEL 

(IE6REES)  (HH/USEC) 

(DEGREES) 

(HH/USEC) 

0  0.71337 

5.35471E-7 

8.71337 

5  8.4039 

.554083 

0.71009 

10  8.59579 

1.12141 

8.70003 

IS  8.44995 

1.70591 

8.48242 

20  0.24794 

2.3209 

8.45478 

25  7.99194 

2.9002 

0.42081 

30  7.48491 

3.70159 

8.57214 

35  7.33048 

4.50927 

8.50489 

40  4.93327 

5.43812 

8.41913 

45  4.49904 

4.54179 

0.29952 

50  4.03515 

7.90947 

0.13244 

55  5.55142 

9.70745 

7.8913 

40  5.04204 

12.2915 

7.52273 

45  4.5902 

14.5943 

4.91422 

70  4.17914 

25.544B 

5.B54B3 

75  3.89977 

45.144 

4.49444 

80  3.77398 

47.7915 

3.8413 

85  3.73184 

01.2412 

3.7398 

90  3.72217 

90.0001 

3.72217 

SLOU  HAVE  SPEEDS  VERSUS  AH6LE 

THETA  PHASE  1  PHASE  VEL 

THETA  GROUP 

S  GROUP  VEL 

(IE6REES)  (HH/USEC) 

(DEGREES) 

(HH/USEC) 

0  2.41549 

8.44539E-4 

2.41549 

5  2.42354 

0.91947 

2.42971 

10  2.4448 

17.5449 

2.44991 

IS  2.48425 

25.4413 

2.73122 

20  2.73414 

33.0443 

2.00481 

25  2.79418 

39.7442 

2.00941 

30  2.84177 

45.7504 

2.97341 

35  2.93409 

51.0419 

3.05328 

40  3.00021 

55.7444 

3.1255 

45  3.08108 

59.0259 

3.10719 

SO  3.14937 

43.2574 

3.2354 

SS  3.20903 

45.0411 

3.24735 

40  3.25407 

44.9417 

3.27024 

45  3.27294 

44.7922 

3.27290 

70  3.23907 

54.9942 

3.35339 

75  3.1050 

39.0444 

3.79770 

•0  2.00531 

37.0404 

3.09111 

•5  2.49095 

53.2410 

3.14411 

90  2.41549 

90. 

2.41549 

A-4 
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APPENDIX  B 

PROGRAM  FOR  CALCULATING  LAMINA  STIFFNESS  COEFFICIENTS 
FROM  PLATE  WAVE  SPEEDS 


The  following  computer  program,  written  in  BASIC,  applies  to  the  calculation 
of  the  four  reduced  lamina  stiffness  coefficients  Qn,  Q22,  Qfi6  and  Q12  when 
more  than  four  experimental  wave  speeds  have  been  recorded.  It  assumes  that  at 
each  angle  interrogated  within  the  lamina  both  the  pseudolongituainal  wave  speed, 
Vp,  and  the  pseudoshear  wave  speed.  Vs,  have  been  measured.  The  program  is 
based  upon  a  least-square-error  fit  of  two  theoretical  curves,  Xp  (e)  and  X5  (e), 
to  the  set  of  experimental  data,  where  X  *  p  V  (e)2.  Enter  the  data  beginning 
in  Step  7000  in  the  order:  angle  of  the  transducer  with  respect  to  the  fibers, 
e  (degrees);  Xp  (GPa)  obtained  at  this  angle;  X$  (GPa)  obtained  at  this  angle. 

The  program  accommodates  three  trial  values  for  each  of  the  four  Qij'S  and 
computes  the  square-error  sum  for  all  81  combinations  of  Qij.  The  output  lists 
the  best  five  combinations  in  order  of  descending  total  square  error;  the 
least-square-error  value,  representing  the  best  fit  of  the  two  curves  simul¬ 
taneously,  is  the  final  set.  Enter  the  trial  Qij 1 s  in  GPa  units  beginning  in 
the  Step  7000  in  the  order:  Q^;  Q22;  066 •  Ql2 -  A  sample  of  the  program  output 
is  presented  in  Table  B-l;  the  numbers  therein  are  representative  of  B/Al.  The 
two  least-square-error  velocity  curves  are  shown  with  the  experimental  data  on 
Figure  B-l. 
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100  REM  PROGRAM  TO  CALCULATE  PLATE  flIJ  PARAMETERS 
110  REM  FROM  SYMMETRICAL  LAMI  AND  SM  SHEAR  HAVE 
120  REM  SPEED  SATA.  INPUT  ANGLE ( DEGREES ),XF(GPA) 

130  REN  XS(SPA).  INPUT  ESTIMATED  VALUES  FOR 

140  REM  Oil,  022,  044,  012  IN  8PA  UNITS. 

990  BIN  IC100), C1(100),C2(100),C3(100),C4(100> 
ffS  REM  N  *  NUMBER  OF  AN8LES  INTERR06ATED 
1000  FOR  N-1  TO  5 
1003  READ  T1(N),X1CN),X2(N) 

1010  T(H)«T1(N)*(3. 14159/180) 

1013  NEXT  M 

1020  FOR  H*1  TO  3 

1023  READ  01 (H), 02(H), 04(H), 03(H) 

1030  NEXT  H 

1031  PRINT  "TRIAL  VALUES  OF  OIJ  IN  GPA  UNITS" 

1032  PRINT  "Oil", "022", "044", "012" 

1033  FOR  H  •  1  TO  3 

1034  PRINT  01(H),02(H), 04(H), 03(H) 

1034  NEXT  H 

1038  M*1 

1040  FOR  1-1  TO  3 
1045  FOR  J-1  TO  3 
1050  FOR  K-1  TO  3 
1055  FOR  L>1  TO  3 

1058  02-0 

1059  REM  N  -  NUMBER  OF  ANGLES  INTERROGATED 
1040  FOR  N-1  TO  5 

1100  A1(N)  «  01 (I)*(C0S(T(N)))**2  ♦  04(L)*(SIN(T(N) ))**2 
1110  A2(N)  -  04(L)*(C0S(T(N) ) )**2  ♦  02(J)*(SIN(T(N)))**2 
1120  A3(N)  «  (03(K)+04(L))  *  (SIN(T(N) )*COS(T(N) ) ) 

1150  X5(N)  •  ( A 1 (N)+A2(N))/2 

1155  X4(N)  «  0.5*(((A1(N)-'A2(N))m2  +  4*A3(N)**2))**0.5 
1140  X3(N)  •  X5(N)  ♦  X4(N) 

1145  X4(N)  -  X5(N)  -  X4(N> 

1180  DKN)  «  ( X3 ( N ) -XI  (N)  )**2  ♦  (X4(N)-X2(N>)**2 
1190  B2-D2*D1(N) 

1200  NEXT  N 
1210  B(M)  «  D2 
1212  CUM)  «  OKI) 

1214  C2(M)  •  02 ( J) 

1214  C3(N)  -  03 ( K ) 

1210  C4(H)  •  04(L) 

1220  N  «  N+1 
1245  NEXT  L 
1230  NEXT  N 
1253  NEXT  J 
1240  NEXT  I 
1400  N-1 

1410  FOR  1-1  TO  00 

1420  N-I  ♦  1 

1.430  FOR  J-K  TO  01 

1440  IF  BCD  >  D(J)  THEN  1500 

1450  T1  «  1(1) 
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14S2  T2  •  CHI) 

1454  T3  -  C2(I) 

1454  T4  •  C3(I) 

1451  TS  «  C4(I) 

1440  B(I)  «  l(J) 

1442  CHI)  >  CHJ) 

1444  C2(I)  >  C2( J> 

1444  C3(I)  >  C3< J) 

1441  C4(I)  >  C4( J) 

1470  l(J)  >  T1 
1472  CHJ)  -  T2 
1474  C2( J)  -  T3 
1474  C3( J)  >  T4 
1470  C4( J)  >  T5 
1500  NEXT  J 

1110  NEXT  I 

1402  PRINT 

1403  PRINT 

1404  PRINT  "REST  FIVE  SETS  OF  Q1J  (CPA)'1 

1405  PRINT 

1404  PRINT  *011", "022" , "844" t  *012", "SQUARE  ERROR" 

1407  PRINT 

1408  FOR  I  >  77  TO  81 

1410  PRINT  CKI),C2(I),C4(1),C3(I),D(I) 

1420  NEXT  1 

1410  REN  ENTER  ANGLE,  XF,  XS 
7000  RATA  00,239.52,58.124 
7100  RATA  20,219.21,41.055 
7200  RATA  45,171.43,72.441 
7300  RATA  40,143.28,71.930 
7400  RATA  90,158.38,40.079 
7450  REN  ENTER  011,022,044,012 
7500  RATA  235,157,54,37.5 
7400  RATA  235.5,157.5,54.5,38 
7700  RATA  234,158,57,38.5 
9999  END 
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TABLE  B-l 

.  REDUCED  STIFFNESS  COEFFICIENTS 

OF  B/A£  from  PLATE  WAVE  DATA 

TRIAL  VALUES  OF  OIJ  IN 

GPA  UNITS 

Oil 

022 

066 

Q12 

235 

157 

56 

37.5 

235.5 

157.5 

56.5 

38 

236 

158 

57 

38.5 

BEST  FIVE 

SETS  OF  OIJ  (SPA) 

on 

022 

066 

012 

SQUARE  ERROR 

236 

157.5 

56.5 

37.5 

90.6936 

235.5 

157.5 

56.5 

38.5 

90.6852 

236 

157.5 

56.5 

38 

90.6342 

235.5 

157.5 

56.5 

37.5 

90.5591 

235.5 

157.5 

56.5 

38 

90.4256 
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